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I. INTRODUCTION 



Exact Expressions for Minor Hysteresis Loops in the Random Field Ising Model on a 

Bethe Lattice at Zero Temperature 

Prabodh ShuklsE 
Physics Department 
North Eastern Hill University 
Shillong-793 022,India. 

We obtain exact expressions for the minor hysteresis loops in the ferromagnetic random field Ising 
model on a Bethe lattice at zero temperature in the case when the driving field is cycled infinitely 
slowly. 
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Hysteresis is a common phenomenon but uncommonly difRcult to treat analytically. Theoretically, hysteresis is 
(y-j ' expected to vanish as the frequency of the driving field goes to zero, or its period goes to infinity. However, in many 
cases hysteresis persists over the longest experimental time scales. For example, there are indications that hysteresis 
would be observed in a permanent magnet even if the applied field were to be cycled so slowly as to take the entire life 
of an experimentalist to complete one loop. It is of practical importance to make a theory for this kind of hysteresis 
which persists over the longest practical time scales, and where the effect of temperature on the area and the shape 
r of the hysteresis loop is small. A reasonable starting point for the theory may be the Glauber dynamics of the Ising 
^ 1 model at temperature T, driven by a field of frequency lo. However, this is difficult to treat analytically. A simpler 
O ' version of this model which appears to be adequate for our purpose is obtained by taking the limit T = first, 
and then the limit w = 0. The zero-temperature, zero-frequency version produces realistic looking hysteresis loops if 
one incorporates a Gaussian distribution of on-site quenched random field in the model. Thus the zero-temperature 
single-spin-flip dynamics of the ferromagnetic random field Ising model (RFIM) was proposed as a model of hysteresis 
and Barkhausen noise by Sethna et. al. Q. Anti- ferromagnetic RFIM is also interesting in the context of relaxation 
without Barkhausen noise 0. The difficulty in the analytical treatment of the models comes from the presence of 
Q the quenched random field. The zero-temperature dynamics of RFIM can not be solved exactly in two or three 
O ' dimensions (so far). We have obtained the major hysteresis loop for the ferromagnetic RFIM in one dimension |3] and 
on a Bethe lattice The Barkhausen noise on the major loop has been analyzed Minor hysteresis loops in the 
ferromagnetic RFIM have been obtained in one dimension 6J. Anti-ferromagnetic RFIM is apparently more difRcult, 
^ , and its analytic solution is limited so far to the major hysteresis loop in one dimension for a rectangular distribution 
' of the quenched field of width 2 A where A < | J|, J being the strength of the nearest neighbor interaction. 

In the following, we present the solution of minor loops in the ferromagnetic RFIM on a Bethe lattice. This is an 
. extension of the work presented in reference , and the completion of a problem which remained open in reference 0] ; 
' therefore the reader is assumed to be familiar with these references. Several important aspects of the analysis and 
simulations of the model which are discussed in references Q and |^ are not repeated here in order to save space, or 
mentioned only briefly for the sake of completeness. It was shown in one dimension that a reversal of the applied 
field by an amount 2 J from anywhere on the major loop brings the system on the opposite half of the major loop. 
' The reversed trajectory merges with the opposite half of the major loop for the portion of the reversed field exceeding 
2J. Thus, the width (along the applied field axis) of a minor loop which touches both halves of the major loop (but 
does not merge with cither of them) is constant and equal to 2 J irrespective of the position of the minor loop. The 
' O ' shape of the minor loop does depend upon its position inside the major loop. This physically interesting result is now 
^ . shown to hold on all Bethe lattices irrespective of the coordination number z of the lattice. 
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II. MAJOR HYSTERESIS LOOP 



In this section, we recall the model briefly, and the solution of the major hysteresis loop obtained earlier The 



5^ ^ RFIM in an external field h^xt is characterized by the Hamiltonian: 
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H = - J ^ SiSj - ^ hiSi — hext ^ Si 
i,j i i 



(1) 



The sum in the first term is restricted to pairs of nearest neighbors on a Bethe lattice of coordination number z. 
The external field hext is cycled from — oo to +00 and back to —00. This takes the system around its major hysteresis 
loop. Spins turn up on the lower half of the loop, and turn down again on the upper half. The applied field changes 
very slowly. Equivalcntly, at each value of the external field, the system is allowed adequate time to attain a relaxed 
state with each spin pointing along the net field at its site. In the relaxed state at hext = h, the probability that an 
arbitrary site i is up is given by, 



p^''^ = E u i^^^'^)]"!! - p*ih)r->n{h) (2) 

n=0 ^ ^ 

Here P*{h) is the conditional probability that a nearest neighbor of site i is up before site i is relaxed , and Pn{h) 
is the probability that the site i with quenched field hi can turn up in applied field /i if n of its nearest neighbors are 
up. The starting state on the lower half of the major loop has all sites down, and the starting state on the upper half 
has all sites up. Thus, on the lower half of the major loop, P*{h) denotes the conditional probability that a nearest 
neighbor of site i turns up before site i. On the upper half, P*{h) denotes the conditional probability that a nearest 
neighbor of site i turns down after site i. We distinguish the two situations by putting a subscript I for the lower half 
and u for the upper half. This gives 



n=0 



Pt{h) = E ( ' n M - Pt{h)Y-'->n{h), 



(3) 



and 



p*u{h) = E n [p:m^ - p*u{h)Y-'->n+i{h) (4) 

n=0 ^ ^ 

We note that pn+i{h) = Pn{h + 2J), and therefore P*{h) = P*{h + 2J). Here, P*{h) is the conditional probabihty 
that given a site which is up on the upper half of the major loop at field h, its nearest neighbor is also up. Similarly 
P;*(/i + 2 J) is the conditional probability that given a site which is down on the lower half of the major loop at field 
h + 2 J, its nearest neighbor is up. 



III. STARTING A MINOR LOOP 



Suppose we are on the lower part of the major loop when the applied field is reversed from h to h' {h' < h) to 
generate the upper half of a minor hysteresis loop. We ask the question, what is the probability that an arbitrary 
site i which was up at h turns down at h' . In order to compute this probability correctly, we must take into account 
the irreversibility of the zcro-tcmpcraturc dynamics. Consider a site i which is down on the lower half of the major 
loop at an applied field h — 6h but turns up at h, where 6h is an arbitrarily small field. Once site i has turned up, it 
may not turn down immediately if the field is rolled back to the value h — 6h. The reason is as follows. When site i 
turns up, it increases the net field on each of its nearest neighbors by an amount 2 J. The increased field may cause 
one or more nearest neighbors of site i to turn up. Suppose Ua nearest neighbors of site i were already up before site 
i turned up, and nt nearest neighbors turn up after site i turns up. Clearly, rih must lie in the range < nt < z — Ua- 
The rib neighbors increase the local field at site i by a finite amount 2n(, J. Therefore an infinitesimal decease in the 
applied field will not cause site i to turn down unless = 0. A site i with m, > will turn down in decreasing applied 
field only after all of the rib nearest neighbors have turned down. When the field is reversed to h' < h, none of the Ua 
neighbors (which turned up before site i turned up) could possibly turn down before site i turns down. This leaves 
the other rift neighbors which turned up after site i. The Ub neighbors turned up because the field on them increased 
by an amoimt 2 J after site i turned up. In decreasing field h' , the rib neighbors will turn down before site i turns 
down. At h' = h — 2 J, all of the rib neighbors would have turned down with the result that the nearest neighbors of 



a site i which are up at /i' = /i — 2 J are precisely those which were up before site i flipped up. These neighbors will 
turn down on the reverse trajectory only after site i turns down. Thus, given an up site i a,t h' = h — 2J on the upper 
half of the minor loop, the conditional probability P^ih — 2 J) that a nearest neighbor of i is up is equal to P*{h), 
where P*{h) is the conditional probability that a nearest neighbor of site i is up at field h given that site i is down 
on the lower half of the major loop at field h. The probability that the site i is up at /i — 2 J is given by the equation. 



Pih-2j) = j2('n) [PtihTil - Pt{h)ryn{h - 2J) 



(5) 



or, using the identity = P*{h- 2 J), 



- 2j) = ^ ( M [p:ih - 2j)]"[i - p:{h - 2j)]^->„(ft - 2j) 

n=0 ^ ^ 



(6) 



It follows from the above equation that the reverse trajectory will meet the upper half of the major loop at 
h' = h — 2 J and merge with it for h' < h — 2 J. Thus, the task of computing the minor hysteresis loop is reduced to 
range h — 2J < h' < h. We return to the question asked at the beginning of this section. What is the probability 
that a site i which is up at h turns down at h' ? This is given by, 

^'e^') = i(n) [Ptih)r[D*{hT-'' [Pn{h)-Pn{h')] (7) 

Here D*{h') is the probability that a nearest neighbor of site i turns down on the reverse trajectory before site i. 
D*{h') is determined by the equation, 

D*{h') = E ( ^ n ^ ) [Ptm^ - Pt{h)Y-'-- [1 -Pn+i{h)] 



n=0 

z-l 



+E ( ) [pt{h)nD*{h')r-'-- \pn+,{h)-pn+,{h')] (8) 

Given a site i which is up at h, the first sum above gives the conditional probability that a nearest neighbor of site 
i is down at h' = h, i.e. at the very start of the reverse trajectory (and hence remains down for h' < h). The second 
sum takes into account the situation that the nearest neighbor in question is up at h, but turns down before site i 
turns down on the return loop. Note that all the nearest neighbors of a site i which turned up after it turned up on 
the lower major loop must turn down before site i turns down on the upper minor loop. 

The magnetization on the reverse trajectory is given by, 

m'{h') = 2\p{h)-q'{h')]-l (9) 

IV. COMPLETING THE MINOR LOOP 

We reverse the field h' to h" {h" > h') to trace the lower half of the return loop. The magnetization on the lower 
half of the return loop may be written as, 

m"{h") = 2\p{h) - q'{h') + p"{h")] - 1 (10) 

where p"{h") is the probability that an arbitrary site i which turned up at h and turned down at h' , turns up again 
at h". 

P"{h") = i2(n) [U*{h"r[D*{hr \Pn{h") -Pn{h')] (11) 



Here U*{h") is the conditional probability that a nearest neighbor of a site i turns up before site i turns up on the 
lower return loop. It is determined by the equation, 

u*{h") = p*{h) - E ( ^ n M [pnh)r[D*{hT-'-''[Pn{h)-pn{h')] 

+ J2('n^)p*('''T{D*{hT-'~''[Pn{h'')-Pn{h')] (12) 

The rationale behind equation (12) is similar to the one behind equation (8). Given that a site i is down at h', the 
first two terms account for the probability that a nearest neighbor of site i is up at h" > h'. Note that the neighbor 
in question must have been up at h in order to be up at h', and if it it is already up at h' then it will remain up on 
the entire lower half of the return loop, i.e. at h" > h'. The third term gives the probability that the neighboring 
site was down at h' , but turned up on the lower return loop before site i turned up. It can be verified that the lower 
return loop meets the lower major loop at h" = h and merges with it for h" > h, as may be expected on account of 
the return point memory. The exact expressions given above have been checked against numerical simulations of the 
model in selected cases resulting in excellent agreement in all cases which were tested. 



V. CONCLUDING REMARKS. 



The method of calculating the minor loop described above may be extended to obtain a series of nested minor 
loops. The key point is that whenever the applied field is reversed, a site i may flip only after all neighbors of site 
i which flipped in the wake of site i (on the immediately preceding sector) have flipped back. The neighbors of site 
i which remained firm after site i flipped previously do not yield before site i has flipped again on the return loop. 
We have obtained above expression for the return loop when the applied field is reversed from h^xt — h on the lower 
major loop to hext = h' {h — 2J < h' < h), and reversed again from h^xt = h' to h^xt = h" {h" < h). When the 
applied field is reversed a third time from h" to h'" [h'" < h"), expressions for the magnetization on the nested return 
loop follow the same structure as the one on the trajectory from h to h' . Qualitatively, the role of P* on the first leg 
{h to h') is taken up by U* on the third leg {h" to h'") of the nested return loop. 

We conclude by comparing the results obtained above with numerical simulations on Bethe lattices of coordination 
number z = 3 and z = 4, and also contrast these results with those obtained in the one dimensional case 0. Let us 
specifically choose a Gaussian distribution of the random field with mean value zero and variance cr^ . One generally 
expects the solution of an Ising model with nearest neighbor interactions on a Bethe lattice (z > 3) to be qualitatively 
different from its solution in one dimension (z = 2), and to be similar to the mean-field solution of an infinite-range 
model. However, these expectations are not born out in the case of hysteresis in RFIM. The mean-field solution jl| 

does not show any hysteresis for a > o'c(oo) = contrast to this, there is hysteresis on Bethe lattices for all 

values of a. Moreover, the behavior on lattices with z = 2 and z = 3 turns out to be qualitatively similar. For 
lattices with z > 4, the magnetization on each half of the hysteresis loop has a jump discontinuity for a < adz)', 
the jump discontinuity is absent on lattices with z < 3 for any finite value of a. Figure 1 shows the major as well 
as a minor hysteresis loop on a Bethe lattice with z = 3, and = 2. The minor loop starts on the lower half of the 
major loop at h = 1.5 and meets the upper half at h = —.5 as may be expected from the theoretical prediction. As 
the magnetization on the lower half of the major loop is a single valued function of the applied field in Figure 1, the 
point where the applied field is reversed determines the minor loop uniquely. The one dimensional case (z = 2) is 
similar [gj. However, the situation is different for z > 4. For z = 4, we have crc(4) = 1.78 approximately. Figure 2 
shows the major loop for z = 4 and a = 1.70 with a jump discontinuity at a critical field he which is slightly higher 
than unity for a = 1.70. There are two values of the magnetization at he both lying on the lower half of the major 
loop. If we reverse the applied field from the value he, we must specify the state of the magnetization from where 
the return is made giving us two possible return loops originating at he. Figure 2 shows two minor loops starting 
at h — .95 (slightly less than he), and h — 1.05 (slightly greater than he) on the lower half of the major loop. Both 
return trajectories touch the upper half of the major loop when the field has been reversed by an amount 2J as 
expected from the theoretical analysis. The overall agreement between the simulations and the theory is also quite 
good indicating that the model considered here is self-averaging ^ IH |^ . 



[1] J P Sethna, K A Dahmen, S Kartha, J A Krumhansl, B W Roberts, and J D Shore, Phys Rev Lett 70, 3347 (1993). 



fig1.ps 




-3-2-10123 

h 

FIG. 1: Major and minor hysteresis loops in RFIM (J = 1) on a Bathe lattice {z = 3) for a Gaussian distribution of the 
random field with mean zero and a = 2. The minor hysteresis loop is obtained by reversing the applied field from h = 1.5 to 
h' = —.5 and back to h = 1.5. Theoretical result is shown by a continuous line, and symbols show the data obtained from 
numerical simulation of the model. 
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FIG. 2: Hysteresis in RFIM (J = 1) on a Bethe lattice {z = 4) for a Gaussian random field of mean zero and a = 1.70. 
Discontinuities in the major loop vanish above crc=1.78. Two minor loops are shown starting on the lower major loop at 
h = .95 and h = 1.05 respectively. As in Figure 1, theoretical result is shown by a continuous line, and symbols show the data 
obtained from numerical simulation of the model. Note that the minor loops touch the upper major loop when the applied 
field has been reversed by an amount 2J. 



